Astronomy & Astrophysics manuscript no. pn'matter'rev 


February 5, 2008 


(DOI: will be inserted by hand later) 





The structure of Planetary Nebulae: theory vs. practice. 

F. Sabbadin^, M. Turatto^, R. Ragazzoni^, E. Cappellaro'^, and S. Benetti^ 



o 
o 

(N 



> 
m 
oo 
(N 

O 

o 

Oh; 

6 ■ 



^ INAF - Astronomical Observatory of Padua, Vicolo dell'Osservatorio 5, 1-35122 Padua, Italy 

^ INAF - Astrophysical Observatory of Arcetri, Largo E. Fermi 5, 1-50125, Italy 

^ INAF - Astronomical Observatory of Capodimonte, Via Moiariello 11, 1-80131 Naples, Italy 

Received November 21, 2005; accepted January 17, 2010 

Abstract. This paper - the first of a short series dedicated to the long-standing astronomical problem of de- 
projecting the bi-dimensional apparent morphology of a three-dimensional mass of gas - focuses on the density 
distribution in real Planetary Nebulae (and all types of expanding nebulae) . We introduce some basic theoretical 
notions, discuss the observational methodology and develope the accurate procedure for the determination of the 
matter radial profile within the sharp portion of nebula in the plane of the sky identified by the zero-velocity- 
pixel-column (zvpc) of high-resolution spectral images. Moreover, a series of evolutive snapshots is presented, 
combining illustrative examples of model- and true-Planetary Nebulae. Last, the general and specific applications 
of the method (and some caveats) are discussed. 

Key words, planetary nebulae: general- ISM: kinematics and dynamics- ISM: structure 



1. Introduction 

The spatial distribution of gas constitutes the most im- 
portant observational parameter for Planetary Nebulae 
(PNe): when combined with gas dynamics, it supphes 
fundamental information on the mechanisms and physi- 
cal processes driving nebular evolution (mass-loss history, 
wind interaction, ionization, magnetic fields, binarity of 
the central star etc.) and allows us a direct comparison 
of each real object with theoretical evolutionary models, 
hydro-dynamical simulations and photo-ionization codes. 

In spite of this, the whole astronomical literature ded- 
icated to PNe contains only a handful of papers dealing 
with the true spatial structure, which is generally obtained 
by deconvolution of volume emissivity from the appar- 
ent image using Abel's integration equation for spherically 
symmetric systems and "ad-hoc" algorithms for inclined 
axi-symmetrical nebulae (Wilson & AUer 1951, Lucy 1974, 
Soker et al. 1992, Volk & Leahy 1993, Bremer 1995). These 
rare - quite rough and controversial - results stress the 
huge difficulties so far encountered in de-projecting the 
bi-dimensional appearance of a three-dimensional struc- 
ture (AUer 1994). 



* Based on observations made with ESO Telescopes at 
the La Silla Observatories (program ID 65.1-0524), and 
TNG (Telescopio Nazionale Galileo) at La Palma, Canary 
Islands (program AOTlO-10). We extensively apply the photo- 
ionization code CLOUDY, developed at the Institute of 
Astronomy of the Cambridge University (Ferland et al. 1998) . 
Send offprint requests to: F. Sabbadin, sabbadin@pd.astro.it 



Recently, Sabbadin et al. (2005 and references therein) 
have suggested that the dynamical properties of the ion- 
ized gas represent the key for overcoming the stumbling 
block of de-projection in PNe and other types of expand- 
ing nebulae, like Nova and Supernova Remnants, shells 
around Population I Wolf-Rayet stars, nebulae ejected by 
Symbiotic Stars, bubbles surrounding early spectral- type 
Main Sequence stars etc.. 

Historically, the general kinematical rule for the bright 
main shell of PNe goes back to Wilson (1950), based on 
Coude spectra - no image-dcrotator - of 26 targets: the 
high-excitation zones expand more slowly than the low- 
excitation ones, and there is a direct correlation between 
the expansion velocity and the size of the monochromatic 
image. 

Weedman (1968) secured Coude spectra (-1- image- 
derotator) of a sub-set of Wilson's list (10 PNe), covered 
along the apparent major axis. Assuming "a priori" the 
prolate spheroid hypothesis with a/b=1.5, Weedman de- 
rived a typical expansion law of the type T4xp=s(R-Ro), 
where s is the (positive) slope of the expansion velocity 
gradient and Ro the radius at which Vi^xp^O. In most cases 
Ro ~0 or Ro<<Rncb (i.e. the outflow is nearly ballistic, 
Hubble- type). 

The classical papers by Wilson (1950) and Weedman 
(1968) belong to the heroic age of photographic plates. 
Later on, the introduction of linear 2-D detectors at high 
efficiency (CCDs) greatly enhanced flux accuracy in a wide 
spectral range, and the use of echelle spectrographs al- 
lowed the observers to extend the analysis to large, faint 
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objects. Even so, a odd dichotomy marks out the whole 
astronomical literature: high-resolution spectra of PNe are 

used to obtain either the detailed kinematics in a few ions 
(generally, [O III], H I and [N II]) or the "average" nebular 
fluxes. 

In the former case (i. e. kinematics), long-slit observa- 
tions are performed through an interference filter to iso- 
late a single order containing one or two nebular emissions. 
Valuable examples concern the kinematics of: 15 multiple 
shell PNe (Guerrero et al. 1998; single position angle (PA), 
instrumental spectral resolution AV=6.5 km s~^), NGC 
2438 (Corradi ct al. 2000: single PA, AV=4.3 km s^^) 
and MZ 3 (a symbiotic-star nebula), covered at various 
PA by Santander-Garcia et al. (2004, AV=6.0 km s'^) 
and Guerrero et al. (2004, AV=8.0 km s"-'^). We must no- 
tice, however, that the standard observational procedure 
- strongly reducing the useful spectral range by means of 
an interference filter - represents an excessive precaution 
in most cases, since the PN is an emission line object. As 
shown by Benetti et al. (2003) and Sabbadin et al. (2004), 
a PN larger (even much larger) than the order separation 
of the echellograms can be covered with a single exposure 
in the whole spectral range, the only limit (variable from 
one instrument to an other) being the effective superposi- 
tion of spectral images belonging to different orders (and 
not the mere order separation). 

In the latter dichotomous case (average fluxes from 
high-dispersion spectroscopy), line intensities are inte- 
grated over the entire slit (whose length is smaller than 
order separation) and provide "mean" physical conditions 
and ionic and chemical abundances within the whole slice 
of nebula covered by the spectrograph, thus losing the de- 
tailed (i. e. pixel-to-pixel) information on kinematics and 
flux. Representative examples are given by Hyung & AUer 
(1998) and Hyung et al. (2001). 

A courageous effort for overcoming the foregoing im- 
passe has been performed by Gesicki & Zijlstra (2000) 
and Gesicki et al. (2003 and references therein): they 
secured high-resolution spectra of a number of compact 
PNe (including Galactic Bulge objects and nebulae in the 
Sagittarius Galaxy and in the Magellanic Clouds) and 
combined integrated emission profiles of forbidden and 
recombination lines with the photo-ionization code of a 
spherical nebula (Torum model), inferring that accelera- 
tion is a quite common property in the PNe innermost 
layers (i. c. "U" -shaped expansion profile), due to the dy- 
namical contribution by the shocked, hot wind from the 
central star. 

Unfortunately, a detailed comparative analysis 
(Sabbadin et al. 2005) proves that the "U" -shaped 
expansion profile is a spurious, incorrect result caused by 
a combination of unsuited assumptions, spatial resolution 
and diagnostic choice. 

All this confirms that, due to the nebular complexity 
and large stratification of the radiation and kinematics, 
the recovery of the spatio-kinematical structure needs ob- 
servations: 



(a) at adequate "relative" spatial (SS) and spec- 
tral (RR) resolutions (SS=R/Ar, R=apparent radius, 
Ar=secing+guiding; RR='\4xp/AV, AV=instrumental 
spectral resolution), 

(b) in a wide spectral range containing emissions at 
low, medium and high ionization and 

(c) at several slit positions over the nebula, 

to be combined with a straightforward and versatile 
method of analysis. 

To this end, in 2000 we started a multi-PA spectro- 
scopic survey of bright PNe in both hemispheres, ob- 
served with ESQ NTT-I-EMMI (80 echelle orders cover- 
ing the spectral range AA3900-8000 A with AV=5.0 km 
s-i) and Telescopic Nazionale Galileo (TNG)-I-SARG (54 
orders, spectral range AA4600-8000 A, AV=2.63 km s"!). 
The pixel-to-pixel analysis of flux and velocity in a large 
number of emissions gives the bi- dimensional ionic struc- 
ture (tomography) of each nebular slice intercepted by the 
spectrograph slit and a 3-D rendering program, assembling 
all tomographic maps, provides the accurate spatial distri- 
bution of the kinematics, physical conditions (To and Ng) 
and ionic and chemical abundances within the nebula. 

Besides the - in fieri - detailed spatio-kinematical study 
of individual PNe (Ragazzoni et al. 2001, Turatto et al. 
2002, Benetti et al. 2003, Sabbadin et al. 2004, 2005), we 
decided to deepen a general, long-lasting open problem: 
inferring the radial density profile in expanding nebulae. 

This introductory, didactic paper contains the detailed 
procedure for the determination of the gas distribution in 
real PNe and is structured as follows: Sect. 2 gives some 
basic theoretical notions, Sect. 3 describes the practical 
application to high-dispersion spectra, Sect. 4 presents a 
series of evolutionary snapshots for model- and true-PNe, 
Sect. 5 contains the general discussion and Sect. 6 draws 
the conclusions. 

In a future, applicatory paper we will compare the 
density profiles observed in a representative sample of 
targets with the expectations coming from recent hydro- 
dynamical simulations and theoretical models, in order 
to disentangle the evolutive phenomenology and physical 
processes responsible for PNe shape and shaping. 

2. Theory 

Following AUer (1984), Pottasch (1984) and Osterbrock 
(1989), the absolute flux emitted in the line A by the ele- 
mentary volume of a steady-state nebula is given by: 

F(X-, A) = AreA^(X")f(X'", A,iVe, T,)ei, (1) 

where Tg is the electron temperature, A'e the elec- 
tron density, X" the ion involved in the transition, 
f(X'",A,iVe,Te) the emissivity function for the line (X'",A) 
and ei the local filling factor (fractional volume actually 
filled by matter with electron density Ne). 

The emissivity function f{X'^ ,X,Ne,Te) depends on the 
excitation mechanism. Recombination through the cap- 
ture of an electron (followed by a cascade to lower lev- 
els) and coUisional excitation by an electron (followed by 
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spontaneous radiation) are the main excitation processes 
in PNe. ^ 

For recombination lines (n=m+l in Eq. (1)) 



f(X-,Aij,iVe,7;) = ^aefi(Aij), 



(2) 



where a^f / (Ay ) is the effective hne recombination coeffi- 
cient (a weak function of Ng) given by 



aeff(Aij) = B(Aij)aeff(Xr) 

and B(Aij), the branching ratio, by 



B(Aij) 



(3) 



(4) 



with Aifc=probabihties for spontaneous radiative decay. 

The emissivity function for collisionally excited hues 
(low-density case; n=m in Eq. (1)) is 



he 

f(X-,A,Te) = yqcoii(A) 

and qco«(A), the collisional excitation rate, 



(5) 



qcou(A) = 8.63 x m-^—)T-°-^e^ , (6) 

where J7=collision strength of the transition averaged over 
the Maxwell distribution, a;i=statistical weight of the 
lower level of the transition, AE=excitation potential of 
the upper level. 

In general, ae//(A) has a weak, inverse dependence 
on Te, whereas qcoH(A), according to Eq. (6), is a strong, 
direct function of To. 

Eq. (1) - relating the local electron density with an 
observable quantity (i. e. the local flux) - represents the 
starting point for the practical determination of the ra- 
dial matter profile in real-PNe (after de-projection of the 
apparent bi- dimensional morphology, as discussed in the 
next section). 



3. Practice 

Let us consider an optically thick, regularly expanding 
nebula at a distance D, powered by a luminous post-AGB 

star at high temperature. A spectrograph long-slit cen- 
tered on the exciting star intercepts the radial slice of 



^ We overlook other possible excitation mechanisms, namely 
(a) Bowen resonance fluorescence (a photon emitted by one ion 
is absorbed by another and degraded in a cascade of emission 
lines via intermediate levels), (b) starlight and/or nebular con- 
tinuum fluorescence (a strong local UV radiation field pumps 
an ion into an excited state), (c) dielectronic recombinations 
(a core electron is excited by capture of a free electron) and 
(d) charge-exchange reactions (an electron is exchanged during 
the collision of a ion with H or He, the most abundant species). 
Their possible role into a specific application (i. e. abundance 
dichotomy from optical recombination and collisionally excited 
lines) will be introduced in Sect. 5. 



nebula shown in the left-most panel of Fig. 1, whose high- 
resolution spectral image in different ions - also presented 
in Fig. 1 - enhances (a) the large stratification of the ra- 
diation and kinematics, and (b) the blurred appearance 
of H I, He I and He H recombination lines, due to a mix 
of thermal motions, fine-structure and expansion velocity 
gradient across the nebula (more details are in Sabbadin 
et al. 2005 and Sect. 4.1).2 
Moreover, Fig. 1 contains: 

- the central-star-pixel-line (cspl), common to all PA, 

representing the matter projected at the apparent po- 
sition of the central star, whose motion is purely radial; 

- a zero-velocity-pixel-column (zvpc, superimposed to 

the [O I] spectral image at A6300 A) , giving the spatial 
profile of the tangentially moving gas at the systemic 
radial velocity. 

Hereafter we will focus on the zvpc: as clearly shown 

in the left-most panel of Fig. 1, it provides the ionic dis- 
tribution within a sharp, central, well-defined portion of 
the radial slice of nebula selected by the spectrograph slit 
(strictly speaking, the zvpc is marginally affected by the 
adjacent layers of the nebular slice - and, at the same time, 
the zvpc marginally contributes to the emission of the ad- 
jacent layers - since the instrumental spectral resolution 
is not infinite). 

According to Eq. (1), the absolute flux emitted in the 
line A by the nebular volume sampled by a single pixel of 
the zvpc is given by 



E(X™, A)z;vpc = F(X™, A) 

zvpc ^ ^zvpc; 



(7) 



where V^ypc (in cm^) is represented by the parallelepiped 



Vzvpc = S X W X d 



with s, pixel size along the slit, given by 

, , .,,.,/ X D(cm) 

s(cm) = pixel height (arcsec I x , 

^ ^ * ^ ^ 206265 

w, pixel size perpendicular to the slit, by 

D(cm) 



w(cm) = slit width(arcsec) x 



206265 



(8) 



(9) 



(10) 



and d, depth of the zvpc along the radial direction, in 
general by 

_i R(A)(arcsec) D(cm) 

d(cm) = (5V(kms ) x ^ ^ \ , -rr x — ^ (H) 

^ ^ ^ ' Vexp(A)(kms-i) 206265 ^ ' 



and, for 14xp=AxR", by 

(5V(kms~-^) D(cm 
d(cm) = ^ — — X 



A(km s-i arcsec-i) 206265 ' 
(5V(km s~^) being the pixel spectral resolution. 



(12) 



^ The twelve spectral images in Fig. 1 refer to NGC 6741 at 
PA=15° (ESO NTT-FEMMI). 
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Fig. 1. Nebular slice intercepted by a spectrograph slit passing through the central star (left-most panel) and the 
related high-resolution spectral image (flux in arbitrary scale) in 12 ionic species arranged (loft to right) in order of 
increasing ionization potential (IP; approximative range). The cspl is also indicated, as well as a zvpc (superimposed 
to the [O I] A6300 A spectral image) with the corresponding part of real nebula (the rectangle centered on the nebular 
slice). 



The corresponding flux received by an observer is 

I(X™, A)corr 



I(X™, A)obs 
with 

I(X™, A)corr = 



10[c(H/3)xf(A)] 



E(X'°,A), 



47rD(cm)2 ^^"^^ 

where c(H/3) is the logarithmic extinction at H/3 and f(A) 
the extinction coefficient given by Seaton (1979). 

In the ionized nebula being A^(H) ~A^(H~'") and 
TVe c^l.l5A''(H+), Eq. (7) can be written in the form ^ 



Niei cx 



I(X-,A) 



f(X-,A,Ne,Te)^ N(X) 



where icf(X")= -^^ry is the ionization correcting factor. 
Recalling that 



^ icf(X') 



(16) 



i<z 



(z=atomic number of the element X), at last we obtain 
the general expression 



N(X) 



i<z 



f(Xi,A,Ne,Te) 



(17) 



providing N^^i'^ in the zvpc from observable quantities (i. 
e. the pixcl-to-pixel flux in one line for each ionic species 
of a given element). 

In practice, for hydrogen - which has a single ionized 
state - Eq. (17) becomes 

^^^^^ (18) 



^,^^^jI(HI,6563) 



aeff(HI,6563) he ' 

^ In this and the following equations, the proportion- 
ality constant is [47r(206265)^ x A(km s"^ arcsec"^)] 
/ [1.15X(5V(km s~^)x s(arcsoc)x w(arcsec) xD(cm)] 
for 14xp=AxR" and [47r(206265)'^ xKxp(A)(km s"^] / 
[1.15X(5V(km s~^)x_R(A)(arcsec)x s(arcsec)x w(arcsec) 
xD(cm)] for Kxp 7^ AxR". 



or 



(13) 7V2^ioc[ 



I(HI,4861)corr 4861 

aeff(HI,4861) ^ he ^ 



(19) 



or using any other H I Balmer recombination line. 

For helium, the contribution of both ionization stages 
must be considered 



iV^ei oc [C(He+) -h C(He++)] x 



N(H) 
N(He) 



with 



(15) C(He-) = 



5876, 



Q!eff (Hel, 5876) he 



and 



C(He++) = [ 



I(HeII,4686)corr 4686 
aeff (Hell, 4686) ^ 



(20) 



(21) 



(22) 



or using any other pair of He I and He II recombination 
lines. 

For heavier elements, we are forced to select collision- 
ally excited emissions of suitable ionic sequences, due to 
the weakness of recombination lines, large stratification of 
the radiation and incomplete ionic coverage of the echel- 
lograms. As tracers, we adopt 0''+0"'"+0+"'" for the exter- 
nal, low-to-medium ionization layers (IP range from to 
55 eV) and Ar++-|-Ar^+-|-Ar*+ for the internal, medium- 
to-high ionization ones (IP range from 28 to 75 eV). ^ 

From oxygen 

AT^e, oc [0(0°) + C(0+) + C(0++)] x ^ (23) 

Note, however, that ionic species higher than Ar*"*" are ex- 
pected in the innermost layers of a PN powered by a central 
star with T, >100 000 K. In this case, A7005 A and/or A6435 
A of Ar*"*" become poor density diagnostics for the internal, 
highest-excitation regions (A3425 A of [Ne V], IP range 97 to 
126 eV, should be preferable). 
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with 

_ I(OI,6300)corr ^ 6300 

^ ' ^qcoii(OI,6300) he ^' 

+ ^ I(QII,7320)eo.. 7320 
^ ^ ^qcou(OII,7320) he ^ 



(24) 



(25) 



(7320 eorresponding to A7320.121 A, the strongest eom- 
ponent of the [O II] red quartet) and 



f.^Q + + -^ ^ jI(OIII,5007)eorr 5007, 



^qcoii (0111,5007) he 



(26) 



From argon 



N(H) 



Nle, oc [C(Ar++) + C{A.v^+) + C(Ar^+)] x ^ (27) 
with 



C(Ar^+) = [ 
and 

C(Ar4+) 



I(ArIV, 4740)c 



4740, 



qcoii(ArIV,4740) he 



I(ArV, 7005)corr 7005 
^qcoii(ArV, 7005) ~h^^ 



(29) 



(30) 



Besides Eqs. (17) to (30), the profile in the zvpe 
is also given by line intensity ratios of ions in p^ eonfigu- 
ration, hke A6717 A/A6731 A of [S II] for low-ionization 
regions and A4711 A/A4740 A of [Ar IV] for the high- 
ionization ones. However, these diagnostics are generally 
weak in PNe and 7Ve[S II] and iVe[Ar IV] can be derived 
only at (or close to) the corresponding intensity peak. 

Summing up: in the most favourable case (i. e. long- 
slit, wide spectral range, absolute flux calibrated eehel- 
lograms of a PN at known distance), we infer both 
Afe(zvpc) and ei(zvpc) by combining Eqs. (17) to (30) with 
A^e[S II] and Afe[Ar IV] (since (a) chemical abundances 
come from line fluxes integrated over the whole spatial 
and spectral profile, Alexander & Balick 1997, Perinotto 
et al. 1998, and (b) the Tc (zvpc) distribution comes from 
diagnostics of ions in p^ and p'' configurations, like A6584 
A/A5755 A of [N II] for low-ionization regions and A5007 
A/A4363 A of [O III] for the high- ionization ones, Turatto 
et al. 2002, Benetti et al. 2003). 

On the other hand, in the most unfavourable case (i. e. 
flux un-calibrated spectra of a PN at un-known distance) 
Eqs. (17) to (30) provide the "relative" N^e]''^ profile in 
the zvpc. ^ 

In order 



^ Moreover, (a) the combination of Eqs. (15) and (17) gives 
the radial ionization structure of each element (more details 
are in Sect. 4.1) and (b) the accurate He/H, O/H and Ar/H 
radial chemical profiles can be obtained by assembling Eqs. 
(18), (19), (20), (23) and (27); we will deepen this topic in a 
future, dedicated paper. 



- to check the general validity and applicability of our 
Nf. reconstruction method for different physical and evo- 
lutionary conditions of the ionized gas and 

- to tackle the various "practical" problems connected 
with the observations of a true-PN (e. g. finite spatial 
and spectral resolutions and incomplete ionic coverage of 
the echellograms, thermal motions, turbulence and fine- 
structure of recombination lines), 

in the next section we 

a) will expose an arbitrary and imaginary model neb- 
ula to a series of UV stellar fluxes and 

b) will analyse in detail the zvpe-flux profile in a few, 
representative true-PNe. 

4. Theory vs. practice 

According to the most recent theoretical evolutionary 

models and radiation- hydrodynamics simulations, a PN 
is the result of the synergistic effects of ionization and 
fast wind on the gas ejected during the superwind phase 
of an AGB-star (Blocker & Schonberner 1990; Vassiliadis 
& Wood 1994; Blocker 1995; Marigo et al. 2001; Perinotto 
et al. 2004b). 

By combining the large variation of post-AGB stel- 
lar characteristics (temperature and luminosity) with the 
gradual nebular dilution due to expansion, we can sketch 
the evolution of a "normal" PN as follows: in the early 
phases the nebula is optically thick to the UV stellar radi- 
ation, then the gas gradually becomes optically thin until 
the quick luminosity decline of the hot star at the end of 
the hydrogcn-shcU nuclear burning, leading to a recom- 
bination of the external nebular layers; later on, the PN 
enters in a final re-ionization phase, because of the slowing 
down in the stellar luminosity decline. 

The recombination and re-ionization phases are absent 
in a low-mass PN ionized by a low-mass, slowly evolv- 
ing central star, whereas massive PNe powered by a high- 
mass, fast evolving post-AGB star never become optically 
thin to the UV stellar radiation. 

A further complication is represented by a possible PN 
"rejuvenation" : the hydrogen-deficient post-AGB star can 
suffer a late thermal pulse during the motion towards the 
white-dwarf region, mixing and burning hydrogen on a 
convcctive turn-over time scale and producing a "born- 
again" PN (Blocker 1995, 2001; Herwig et al. 1999). 

Of course, all this represents an over-simplification, 
each real PN being an extremely complex and inhomo- 
geneous structure: it can be optically thick at some radial 
directions (e.g. along and close to the dense equatorial 
regions) and thin at other directions (e. g. along the low- 
density polar caps). 

For illustrative purposes, we will expose a single, "stan- 
dard" radial matter profile (representing the gas distribu- 
tion of an ideal, average, spherically symmetric nebula) to 
a series of UV stellar fluxes mimicing the evolution of a 
post-AGB star. The resulting ionization structure of each 
"model-PN" snapshot, provided by the photo-ionization 
code CLOUDY (Ferland et al. 1998), will be analysed in 
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Table 1. Input parameters for the adopted, spherically 
symmetric model nebula (CLOUDY) 



Gas density profile 



The same, standard and arbitrary 
double-peak+halo distribution at 
all the four evolutionary steps 
(see Figs. 2, 4, 5 and 6; cfr. Sect. 4) 



Chemical abundances Hc=-0.958, C=-3.30, N=-3.75, 



log[N(X)/N(H)] 



Dust 

Local filling factor 



O=-3.30, F=-6.522, Ne=-3.90, 
Na=-5.52, Mg=-5.30, Al=-6.568, 

Si=-5.00, P=-6.70, S=-5.10, 
Cl=-6.66, Ar=-5.60, K=-7.11, 
Ca=-6.64, Fe=-6.30, 
other elements= CLOUDY default 

CLOUDY default 

1.0 



detail and coupled with the zvpc-flux profiles observed in 
a true-PN at a comparable evolutionary phase. 

The input parameters for the model nebula arc sum- 
marized in Table 1, whereas Table 2 contains the in- 
tegrated flux in the main optical emissions (relative to 
I(H/3)=100), the absolute H/3 flux and the resulting ion- 
ized mass at four representative evolutionary steps. 

Let us start with a "quite young" PN close to the thick- 
thin transition, ionized by a luminous post-AGB star of 
0.605 Mq at a moderate temperature (T*=50 000 K), thus 
postponing to Sect. 4.4 the special case combining proto- 
PNe and low-excitation "born again" PNe. 

4.1. Early PN evolution (Column 4 of Table 2) 

In this first evolutionary snapshot the ionizing engine 
of the adopted, spherically symmetric matter distribu- 
tion (Table 1) is a hydrogen-burning post-AGB star of 
0.605 M0 with T*=50000 K, L*/L0=65OO and tsvK=time 
elapsed from the end of the superwind^SOOO yr (Blocker 
1995 and references therein). 

The results for the model nebula (summarized in Table 
2, Column 4) are graphically shown in Fig. 2, where 

top panel: standard radial matter profile and resulting 

trend for the physical conditions, 
next to top panel: radial ionization structure of the gas, 
next to bottom panel: absolute flux in the main emis- 
sions, 

1 /2 

bottom panel: Nee/ distribution given by Eqs. (17) to 



1 /2 

(30). In our case A^e^/ =-^e 
ei=l (see Table 1 ). 



since we have assumed 



High ionization emissions, like He II and Ar V, are very 
wciak or absent in Fig. 2 (as expected from the "moder- 
ate" star temperature), whereas the enhancement of [O 
I] and [O II] at the extreme nebular edge indicates that 



the model-PN is nearly optically thin to the UV stellar 
radiation. 

To is high (>10 000 K) in the innermost, tenuous lay- 
ers, rapidly decreases to ~8000 K in the densest regions, 
slightly increases further (~9500 K) and definitively falls 
in the outermost, partially neutral strata. 

1 /2 

The bottom panel of Fig. 2 shows the A'eC; pro- 
file obtained from H+, He+-FHe++, 0°-F0+-F0++ and 
Ar++-FAr3++Ar''+ (through Eqs. (17) to (30) and as- 
suming Te=8000 K for all ions, except in the innermost 
nebular regions, where Te(He++)=ro(Ar4+)=13000 K). 
Please note the satisfactory agreement with the input ra- 
dial density distribution (top panel of Fig. 2), confirming 
the general validity (and applicability) of the procedure 
developed in Sect. 3. ^ 

From theory to practice: NGC 6572 (PNG 034.6-hll.8, 
Acker et al. 1992) is a representative example of PN in a 
quite early evolutionary phase. It consists of a bright (log 
F(H/3)=-9.82 erg cm^^ s^^), moderate excitation (excita- 
tion class 5, Hyung ct al. 1994), high-density (logiVc=3.7 
to 4.4, Hyung et al. 1994, Liu ct al. 2004) irrcgidar disk 
(d~5.0") embedded in a faint envelope (10.0"xl6.0") elon- 
gated in PA~0°. 

Wilson (1950) measured Vc^p{[0 I])=16.0 km s'S 
V;xp([0 II])=16.85 km s'^, V;xp([S II])=15.5 km s'^ and 
FexpGN II])=14.75 km s-\ Weedman (1968) gave the gen- 
eral expansion law T4xp(km s~''^)=3.4xR" and Miranda et 
al. (1999, Ha-|-[N II] echellograms at 3 PA; AV=8.0 km 
s~^) evidenced the presence of a coUimated bipolar out- 
flow along and close to the apparent major axis and of an 
equatorial density enhancement of toroidal structure with 
Ve^p(Ra)=U.O km s"^ and V;xp([N II])=18.0 km s"^. 

The exciting star of NGC 6572 has ihb ~13.0, spectral 
type Of/WR (Acker et al. 1992), log(TzH I)~log(TzHe 
11)^4.83 (Phillips 2003 and references therein) and log 
L,/Lq=3.50(±0.20) (for an assumed distance of 1200 pc, 
Hajian et al. 1995, Kawamura & Masson 1996, Miranda 
et al. 1999, and c(H/3)~0.45, Hyung et al. 1994, Liu et 
al. 2004). It is losing mass at a large, although uncer- 
tain, rate (1.0x10^° Mq yr^^, 6.3x10"^ Mq yr^^ and 
3.0x10"^ Mq yr^^ according to Cerruti-Sola & Perinotto 
1985, Hutsemekers & Surdej 1989 and Modigliani et al. 
1993, respectively) and a terminal wind velocity of 1800 
km s~^. 

The closeness of H I and He II Zanstra temperature 
of the star and the enhancement of low ionization species 
([O I], [O II] etc.) at the nebular edge suggest that NGC 
6572 is still optically thick to the UV stellar radiation. 

We have observed the nebula at six equally spaced 
PA with the high-resolution cross dispersed echelle spec- 
trograph SARG (Gratton et al. 2001) mounted at TNG 
(Telescopio Nazionale Galileo), under non-photometric 



^ Although the adopted matter distribution of the model 
nebula is open to criticisms (it is arbitrary, schematic, static 
etc.), we underline that the choice of more sophisticated and/or 
time dependent density profiles does not modify the general re- 
sults illustrated here and in the following sub-sections. 
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Table 2. Four-step model-PN evolution. Integrated intensity of representative nebular emissions (relative to 
I(H/3)=100), absolute H/3 flux and ionized mass for the adopted, spherically symmetric nebula (Table 1) powered 
by a hydrogen-burning post-AGB star of 0.605 Mq. 



A (A) Ion 1(A) 



Sect. 4.4 Sect. 4.1 Sect. 4.2 Sect. 4.3 

proto-PN + early-PN intermediate-PN late-PN 
"born-again" low-exc.-PN 

(very thick) (nearly thin) (thin) (recombining) 

T*=25000 K T,=50000 K T,=100000 K T*=150000 K 

L*/L0=65OO L*/Lo=6500 L*/L0=65OO L*/L0=1OOO 



3425 
3726+29 
3968 
4686 
4740 
4861 
5007 
5876 
6300 
6563 
6584 
6717+31 
7005 
7135 
7320 



[Ne V] 
[O II] 
[Ne III] 
He II 
[Ar IV] 

H I 
[O III] 
He I 
[OI] 

H I 
[N II] 
[S II] 
[Ar V] 
[Ar HI] 
[O 11] 



0.0 
133.3 
0.02 
0.0 
0.0 
100.0 
5.5 
3.11 
1.3 
297.2 
265.5 
43.1 
0.0 
1.6 
0.8 



0.0 

231.2 
15.2 
0.25 
0.22 
100.0 
490.0 
15.5 
2.5 
291.7 
169.0 
24.4 
0.0 
17.4 
2.2 



4.1 
73.8 

49.8 
19.0 

3.4 
100.0 
1510.0 
12.6 

0.1 
285.7 
37.3 
12.4 
0.77 
21.0 

1.0 



14.7 
587.5 
64.9 
42.2 
2.9 
100.0 
1312.2 
10.5 
26.3 
288.4 
434.5 
107.3 
1.05 
30.7 
9.8 



log E(H/3) (erg s'^) 



34.549 



35.057 



35.019 



34.321 



Mi, 



(Mo) 



0.10 



0.38 



0.40 



0.07^*' 

(*) excluding the ionized mass 
of the recombining halo 



sky conditions and seeing ranging between 0.50" and 
0.70". The spectrograph slit (0.40" wide and 26.7" long) 

was centered on the exciting star. The echellograms (ex- 
posure time 360s) cover the spectral range AA4600-8000 A 
with instrumental resolution AV~ 114000 (2.63 km s~^), 
pixel spectral scale (5V=1.20 km s~^ and pixel spatial scale 
s=0.167 arcsec. The spectra were reprocessed according to 
the straightforward procedure described by Turatto et al. 
(2002), including bias, zero order flat field and distortion 
corrections and wavelength and relative flux ^ calibrations. 

Let us select the spectrum at PA=102°, i. e. close to 
the apparent minor axis of NGC 6572, since, according to 



Only "relative" fluxes axe considered here, since the obser- 
vations have been secured in non-photometric sky conditions. 
Note, however, that an absolute calibration can bo obtained by 
means of the H/3, Ha and [O HI] HST imagery of NGC 6572; 
in fact, the absolute flux within the 0.40" x 26.7" rectangle of 
HST image intercepted by the spectrograph slit matches the 
integrated flux of the corresponding spectral emission. 



the de-projection criteria introduced by Sabbadin et al. 
(2005), Rzvpc (apparent minor axis)~ Rcspi- 

The cspl vs. zvpc relation for NGC 6572 at PA=102° 
- based on ten ionic species, i.e. H+, Hc+, O", O^, 0++, 
N+, S"*", S"*"^, Ar++ and Ar+'^ - gives the expansion law 
ycxp(km s~'^)=8.0(±I) xR", whose slope is much larger 
than Weedman's (1968) report. * This implies that the 
dynamical age of the nebula is about 750 yr and the true 
age 1200-1400 yr. 

Due to the nebular compactness and low expansion 
velocity (SS~RR~5, see Sect. 1), the spectral images 
must be carefully de-blurred. To this end we apply the 
Richardson-Lucy algorithm (Richardson 1972, Lucy 1974) 
with point-spread function given by a bi-dimensional 
Gaussian profile characterized by Wgeeing and W(vel)tot5 
where: 



® Weedman obtained V'exp(km s '^)=3.4xR" assuming 
Rzvpc (apparent major axis) ~ Rcspi- 
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Fig. 2. Early modcl-PN evolution (Column 4 of Table 2). General properties of the standard, spherically symmetric 
model nebula (Table 1) powered by a central star with T, =50 000 K and L*/L0=65OO. Top panel: adopted radial 
matter profile and resulting trend for the physical conditions. Next to top panel: radial ionization structure of the gas. 
Next to bottom panel: absolute flux in the main emissions. Bottom panel: Ngc]^'^ distribution given by Eqs. (17) to 
(30), assuming To=8GGG K for all ions, except in the innermost nebular regions, where Tc(He++)=rc(Ar^+)=13 000 

1/2 

K. NqC^ can be directly compared with the electron density profile shown in the top panel, since we have adopted 
e; (model nebula)=l. 

Wsocing - FWHM of the instrumental profile along the Wthormai = thermal motions = 21. 6xlO~^xTo"'^x m~°-^ 
spatial axis of the echellograms - corresponds to 0.52" km s~^ (m=atomic weight of the element); we assume 
(3.11 pixels), Te=11000 K, 

W(vel)tot - FWHM of the broadening function along the Wturb. = random, small-scale motions = 5.0 km s~"^ 
velocity axis of the echellograms - is given by (adopted value), 

Wfine-s. = fine-structure of recombination lines = 7.7 km 
W(vel)tot = [W^ARG + W^hermai + ^tnrb. + WLe-s.]^^^(31) s'^ for H/3 and Ha and 5.0 km s"! for A4686 A of He 

II (after suppression of the 20 km s~^ blue-shifted tail; 
with for details, see Sabbadin et al. 2005) and A5876 A of 

He I. 

WsARG = instrumental resolution = AV =2.63 km s , 
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Fig. 3. Early evolution of a true-PN. General properties in the zvpc of NGC 6572 at PA=102° (apparent minor axis). 
Top panel: flux distribution of the main nebular emissions across the whole nebula (in arbitrary scale and corrected for 
c(H/3)=0.45). Middle panel: re([N II]) (circles) and Ne{[S II]) (squares) at the corresponding line peak and ionization 

1 /2 

structure of oxygen (from Eq. (32) and assuming Te= 11000 K all across the nebula). Bottom panel: relative Ne€i 
profile from A6563 A (H+), A5876 A (He+), A63GG A (O") + A732G A (0+) + A5007 A (0++) and A7135 A (Ar++) 
+ A4740 A (Ar^+) (through Eqs. (17) to (30)). To be noticed: the zvpc crosses the entire true-PN and provides the 
gas characteristics at two opposite radial directions in the plane of the sky. The central star position is at r=0.0 arcsec 
and the slit orientation is shown in the top panel. 



W(vel)tot amounts to ~23 km s ^ (~19 pixels) for H/3 
and Ha, ~13 km s'^ (~11 pixels) for A5876 A of He I 
and A4686 A of He II, and ~7 to 8 km s'^ (6 to 7 pixels) 
for the forbidden lines of heavier elements (N, O, S, Ar 
etc.). The detailed results for the zvpc of NGC 6572 at 
PA=102° are illustrated in Fig. 3, whose top panel con- 
tains the flux distribution across the whole nebula ^ (in 

® According to Fig. 1, the zvpc of a true-PN crosses the entire 
object and provides two independent density profiles at oppo- 



arbitrary scale and corrected for c(H/?)=0.45) of the main 
nebular emissions. Note the absence of the highest ioniza- 



site directions (in Fig. 3 the central star position is at r=0.0 
arcsec and the slit orientation is indicated in the top panel). 
Of course, the two zvpc gas distributions should coincide for 
a spherically symmetric nebula (as in the case of the standard 
model-PN adopted in this paper; that's why Fig. 2 and the 
next figures concerning the model-PN show the radial density 
profile at a single direction). 
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tion lines: A6435 A and A7005 A of [Ar V] and A6560 A 
of He II are below the detection limit, whereas A4686 A 
of He II lies at the extreme blue edge of the frame and is 
too noisy. 

The middle panel of Fig. 3 shows: 

(a) Te([N II]) (circles) and N^{[S II]) (squares) at the 
corresponding line peak. Both diagnostics are uncertain, 
since the [N II] auroral line is quite faint and the A6717 
A/A6731 A intensity ratio close to the high-density limit, 
where collisional de-excitations dominate and the diagnos- 
tic is a weak function of N^; 

(b) the radial ionization structure of oxygen. 
According to Eqs. (15) and (17), it is given by 

I(X',A)oo„ 

^(X') ^ 1 ^ f(X-,A,No,Te) , . 

N(X) icf(Xi) V - 

with i=0, 1 and 2 (moreover, we take To= 11 000 K in the 
whole nebula). Although Eq. (32) is stricty valid in the 
portion of zvpc containing the emission of all the three 
ionic species, we can assume ^^^j =0 for I(O°,6300)=0, 
I(O+,7320)>0 and I(O++,5007)>0. i° 

1 /2 

Fig. 3 (bottom panel) shows the relative N^Ci profile 
from H+, He+, 0°-(-0+-hO++ and Ar++-hAr3+ (through 
Eqs. (17) to (30) and assuming 11 000 K all across the 
nebula). As expected, the ionic sequence of argon under- 
estimates Nosy^ in the external, low ionization regions. 

1/2 • 

Moreover, A^eC/ given by hydrogen is weakened by the 
complex de-blurring procedure and oxygen slightly under- 

1/2 

estimates A^eC; in the innermost, highest ionization lay- 
ers (the same is valid for helium and argon, due to the 
absence of He II and Ar V emissions in the spectrum) . 

Despite these minor discrepancies, the bottom panel of 
Fig. 3 accurately reconstructs the broad and asymmetrical 

1/2 

"relative" Nce^ profile along the entire true minor axis of 
NGC 6572; when combined with Ne{[S II]) (middle panel), 
it provides Ae(top)~21 000 cm'^ (for e/=l). 

The matter distribution at the other PA confirms the 
inhomogeneous and elongated structure of NGC 6572, 
consisting of a dense equatorial torus and two extended 
and fast polar caps at lower density, embedded in a ten- 
uous, ellipsoidal, low-excitation cocoon. The nebula is 
optically thick to the stellar radiation and the ionized 
mass, given by the observed A'e spatial distribution, the 
H/3 flux, and the radio flux (Aller 1984, Pottasch 1984, 
Osterbrock 1989, and Turatto et al. 2002), amounts to 
Mion^0.08(±0.02) M0. 

We end recalling that 

(a) the overall spatial structure of NGC 6572 is avail- 
able at http:/ /web. pd.astro.it/sabbadin, showing the 
multi-color appearance and opaque reconstruction in [O 

Please note that the radial ionization structure given by 
Eq. (32) is independent on the nebular spatio-kinematics (i. e. 
expansion velocity field and matter distribution), since both 
the numerator and denumerator refer to the same pixel of the 
zvpc. 



HI] and [N II] (at high, medium and low flux cuts) for a 
rotation around the East- West axis, 

(b) further examples of "quite young" PNe ionized by a 
luminous star at moderate temperature, close to the thick- 
thin transition are NGC 6210, NGC 6567, NGC 6803, 
NGC 6826 and NGC 6891. 

4.2. Intermediate PN evolution (Column 5 of Table 2) 

In this second evolutionary snapshot the usual, standard, 
double-peak-|-halo radial matter profile of the spherically 
symmetric model nebula (Table 1) is ionized by a 0.605 
Mq hydrogen-burning post-AGB star with T» = 100 000 K, 
L,/L0=65OO and tsw =^4900 yr (Blocker 1995). 

Column 5 of Table 2 gives the overall properties of 
the model nebula, whereas the detailed radial character- 
istics (adopted density profile, physical conditions, ion- 
ization structure and electron density profile obtained by 
Eqs. (17) to (30)) are shown in Fig. 4. The strength of He 
II and [Ar V] emissions, the absence of [O I] and weakness 
of [O II] indicate that this high-excitation model-PN is 
decidedly thin to the UV steUar flux. exceeds 20 000 K 
in the innermost, low-density layers, gradually decreases 
to ~10 000 K in the main regions and increases to ~11 000 
K further. 

Going to the bottom panel of Fig. 4, the Ha flux, com- 
bined with the assumption Tc(H+) = 10 000 K all across the 
nebula, gives a satisfactory N^el^^ profile; it can be fur- 
ther improved adopting distinct values for the inner 
and the outer layers. 

The same consideration applies to helium lines (for 
rc(He++)=15 000 K and To(He+)=10 000 K). 

The ionic sequence of oxygen essentially reduces to 
(with a minor contribution of 0+) and closely re- 
produces the external, medium excitation regions of the 
model nebula, but strongly underestimates the internal, 
high-excitation ones (adopting rc(O+)=rc(O++) = 10 000 
K). 

Argon provides a satisfactory matter distribution 
across the whole nebula, except in the innermost regions, 
where Ar^"*" and even higher ionization species prevail (in 
Fig. 4, bottom panel, we have used ro(Ar'*+)=15 000 K, 
re(Ar3+)=12 000 K and Tc(Ar++) = 10 000 K). 

From theory to practice: NGC 7009 at PA=169° (ap- 
parent minor axis) is the representative example of an 
optically thin PN powered by a luminous post-AGB star 
at high temperature, as shown in the extensive study by 
Sabbadin et al. (2004), based on the absolute flux cali- 
brated ESQ NTT-I-EMMI echellograms at 12 PA. We re- 
fer the reader to this paper for a practical application 
of the model-nebula discussed in this section (as well for 
a detailed reconstruction of the intriguing 3-D ionization 
structure of the Saturn Nebula) . 

More examples of optically thin, intermediate evolu- 
tion PNe excited by a luminous star at high temperature 
are: NGC 2022, NGC 4361, NGC 6058, NGC 6804, NGC 
7094 and NGC 7662. 
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Fig. 4. Intermediate model-PN evolution (Column 5 of Table 2). Radial characteristics of the usual, standard, spheri- 
cally symmetric model nebula (Table 1) ionized by a 0.605 Mq hydrogen-burning post-AGB star with T*=100000 K 
and L«/L0=65OO. Same symbology as Fig. 2. 



We end noting that low-mass PNe ionized by a low- 
mass, slowly evolving central star become optically thin- 
ner and at higher excitation in the late evolution (thus, 
the considerations developed in this section remain valid) , 
whereas "normal" and/or massive PNe powered by a 
"normal" and/or massive star turn to a recombination- 
reionization phase (as discussed in the next section). 

4.3. Late PN evolution (Column 6 of Table 2) 

According to Blocker (1995), at time tsw —7400 yr 
the 0.605 M0 hydrogen-burning post-AGB star has 
T* = 150 000 K and L*/Lq=1000, i. e. it is exhausting 
the hydrogen-shell nuclear burning and quickly fading 



along the white dwarf cooling sequence. The decreasing 
UV stellar flux is unable to fully ionize the nebular gas 
and the outer regions recombine, thus producing a faint 
halo (Tylenda 1986, Phillips 2000, Turatto et al. 2002, 
Sabbadin et al. 2005). 

Although the ionization and thermal structure of this 
model-PN are out of equilibrium, the static assumption 
remains - to a first approximation - valid, since the recom- 
bination time, trec=l/(aB X A'c) (Q!B=cffective recombina- 
tion coefficient), is short: less than a century for hydrogen 
and even shorter for higher ionization species. Of course, 
the steady-state photo-ionization code CLOUDY cannot 
take into account the presence of a recombining halo; this 
it shown in Fig. 5 - containing the radial characteristics 
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Fig. 5. Late modcl-PN evolution (Column 6 of Table 2). General properties of the usual, standard, spherically sym- 
metric model nebula (Table 1) ionized by a star with T* = 150 000 K and L»/Lo = 1000. The presence of a recombining 
halo is ignored by the steady-state photo-ionization code CLOUDY. Same symbols as Figs. 2 and 4. 



of the usual, standard model nebula (Table 1) - and in 
Column 6 of Table 2 - giving the overall intensity in the 
main lines, the absolute H/3 flux and the ionized mass -. 

Note (a) the very-high ionization degree of the inner- 
most regions, (b) the strength of [O I] and [O II] emissions 
at the ionized edge and (c) the temperature radial profile: 
Tg exceeds 30000 K in the innermost, low-density layers, 
rapidly drops to ~11500 K in the densest regions and 
gradually decreases further. 

In Fig. 5 (bottom panel) the line-fluxes of both 
hydrogen (for Te(H+)=11500 K) and helium (for 



Points (a) and (c) explain the central hollow normally ob- 
served in the Ha, [O III] and [N II] imagery of evolved PNe, 
like NGC 6720 (Ring Nebula) and NGC 7293 (Helix Nebula). 



Te(He++)=18000 K and Te(He+)=ll 500 K) closely 

1/2 

reproduce the Nc^i model profile, whereas ar- 
gon (for Te(Ar4+)=20000 K, Te(Ar3+)=15 000 K and 
Te(Ar++)=ll 500 K) underestimates the density of the 
innermost layers (where Ar^+ and even higher ionization 
species dominate) and oxygen (for Te(0"'"+)=11 500 K, 
Te(O+)=9 500 K and Te(O°)=8000 K) only agrees in the 
external, medium-to-low excitation regions of the model 
nebula. 

Going to practice, NGC 6565 (Turatto et al. 2002) and 
NGC 6741 (Sabbadin et al. 2005) are representative exam- 
ples of recombining PNe; both are high-mass objects that 
never became optically thin to the UV radiation of the 
massive, fast evolving central star. The reader is referred 
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to these papers - based on absolute flux calibrated ESO 
NTT+EMMI echellograms at 6 (NGC 6565) and 9 (NGC 
6741) PA - for a practical application of the model-PN 
developed in this section. 

Following the selection criteria introduced by Sabbadin 
et al. (2005), further rccombining PNc arc: NGC 2440, 
NGC 2818, NGC 6302, NGC 6445, NGC 6620, NGC 6886, 
NGC 6894, NGC 7027, IC 4406 and Hu 1-2. 

4.4. A step backward for a special model killing two 
birds with one stone: proto-PNe and 
low-excitation "born-again" PNe ( Column 3 of 
Table 2) 

According to the evolutionary tracks by Blocker (1995), 
the combination T, =25 000 K and L,/L0=65OO identi- 
fies a 0.605 Mq hydrogen-burning post-AGB star at time 
tsw —1600 yr, corresponding to the proto-PN phase: the 
newly born, dense, compact, very-low excitation nebula is 
only partially ionized by the (still modest) UV flux of the 
central star. Typical examples are He 2-47, He 2-131, He 
2-138, Tc 1 and M 1-46. 

On the other hand, the same stellar characteristics 
are found in very-low excitation objects, like NGC 40 
and BD+30°3639, powered by a low temperature and 
high luminosity star of late Wolf-Rayet spectral type and 
widely regarded as "born-again" PNe (Blocker 1995, 2001; 
Herwig ct al. 1999). 

Thus, our standard model nebula (Table 1) exposed to 
the radiation field of a post-AGB star with T* =25 000 K 
and L*/Lo=6500 refers to both the proto-PN phase and 
the low-excitation born-again PN scenario. 

The general properties of the model nebula, summa- 
rized in Tabic 2 (Column 3) and graphically shown in Fig. 
6, can be synthesized as follows: 

- it is optically thick (very thick!), 

- at a very-low excitation degree (note, in particular, the 
weakness of A5007 A ([O III]) and absence of He H, Ar 
IV and Ar V emissions), 

- with a moderate electron temperature: only in the in- 
nermost layers Tq exceeds 8000 K, decreasing to ~6200 
K in the densest regions (with an external bump up to 
7900 K) and declining further. 

1/2 

Concerning the A^oe^ ' profile reconstructed from line 
fluxes through Eqs. (17) to (30) (bottom panel of Fig. 

^'^ Note, however, that collisional de-excitation of forbidden 
lines is important for a typical, very-high density (A'e > 3x10"* 
cm~^) proto-PN. The collisional de-excitation rate is 

gde-ex(A) = 8.63 X 10-^(— )Te-°■^ (33) 

thus, in the high-density case, Eq. (6) becomes 

q,ff{X) = 8.63 X 10-''( — )T-°-5 x (e^ - 1). (34) 



6), a satisfactory agreement is provided by H+ (assum- 
ing Te(H+)=6200 K), whereas He+ (Te(He+)=6200 K) 
and Ar++ (To(Ar++)=6500 K) strongly undcrstimatc the 
density of the main to external layers at low-excitation. 
The ionic sequence of oxygen, 0°-|-0+-|-0++, closely re- 
produces the matter distribution in the inner to main neb- 
ular regions, but presents a large discrepancy in the out- 
ermost parts (for (Te(O++)=6500 K, Te(O+)=6200 K and 
To(O")=6000 K). This curious bump is mainly ascribable 
to a combination of 

(a) low excitation degree of the gas (0+ is the dominant 
species in a large fraction of the ionized nebula), 

(b) steep dependence of qcoii{0 1 1, 7 320) on Tg (much 
steeper than for the other forbidden lines of the sam- 
pie). 

Thus, our assumption Te(O+)=6200 K, constant all across 
the O"^ zone, fails in the external regions, where Tg 
raises up to 7900 K (top panel of Fig. 6): it overes- 
timates N(0+) (and the corresponding A^o(0"'") contri- 
bution) by a factor of two, being [(j'coij(7320, Te=7900 
K)]/ [9coz;(7320,Te=6200 K)]~4.3 (whereas for the other 
forbidden lines of the sample [qcoii{\ Tq=7900 K)]/ 
[qeoii{X,T^=6200 K)]=1.5 to 2.3). 

Passing to practice, an illustrative example of very-low 
excitation PN ionized by a luminous, low temperature star 
is BD-|-30°3639, covered at six equally spaced PA with 
TNG-KSARG. 13 

BD-|-30°3639 (Campbell's hydrogen envelope star, 
PNG 064.7-^05.0 Acker et al. 1992) is a bright (log 
F(H/3)=-10.03 erg cm-^ s'^), high density {\ogNe=3.8 
to 4.3, Barker 1978, Kingsburgh & Barlow 1992, Allcr 
& Hyung 1995, Bernard-Salas et al. 2003), roundish, 
sharp ring (5.0" x 7.0") surrounded by a massive en- 
velope of neutral gas, at a distance of 1200-1500 pc 
(Kawamura & Masson 1996, Li et al. 2002). It is pow- 
ered by a WC9 star with my=12.5, T* ~35 000 K and 
L,/L0 ~4300 (Pottasch ct al. 1978, Koppcn & Tarafdar 
1978, Leuenhagen et al. 1996), which is losing mass at a 
"moderate" wind velocity (~700 km s"^) and a high - 
although uncertain - rate (dc Frcitas-Pacheco et al. 1993 
and Leuenhagen et al. 1996 reported 6.7x10"^ Mq yr"^ 
and 1.3x10"^ Mq yr"^ assuming a nebular distance of 
2.0 Kpc and 2.68 Kpc, respectively). 

Following Bryce & Mellema (1999), BD+30°3639 
presents a peculiar kinematical structure among PNe, the 
internal [O III] layers expanding faster than the exter- 
nal [N II] ones (Kxp([0 III])=35.5 km s"! vs. 14xp([N 
II])=28.0 km s-i). HST STIS observations in the C II] 
line at A2327 A by Li et al. (2002) give yexp(C II])=36.25 
km s~^. 



" The spectra of BD+30°3639 and NGC 6572 have been 
obtained in the same (non-photometric) night and with the 
same instrumental configuration. Thus, the observative infor- 
mation presented in Sect. 4.1 also apply to the echellograms of 
BD-|-30°3639. 
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Fig. 6. Model proto-PNe + low-excitation "born-again" PNc (Colmnn 3 of Table 2) . Radial characteristics of the 
usual, standard, spherically symmetric model nebula (Table 1) for an ionizing star with T*=25000 K and L*/L0=65OO. 
Same symbols as Figs. 2, 4 and 5. 



All this is confirmed and deepened by our 
TNG+SARG echellograms. In general, they indicate 
that: 

(a) the overall excitation degree of the gas is extremely 
low: mean-ionization emissions (c. g. O III, He I and Ar 
III) are weak and high-ionization lines (He II, Ar IV and 
Ar V) totally absent, 

(b) mean-ionization species move faster than the low- 
ionization ones, 

(c) at each PA the range in Rzvpc is quite small, sug- 
gesting a sharp radial density profile, 

(d) faint, fast-expanding knots and tails {Ve^p up to 
100 km s~^) are present in the strongest emissions (as 
already noticed by Bryce & Mellema 1999), 



(e) the spectral images are characterized by a series of 
local, small-scale irregularities and deformations. 

In detail, the cspl vs. zvpc relation at PA=10° (ap- 
parent minor axis of BD-|-30°3639), based on nine ionic 
species, enhances the presence of two distinct velocity 
fields: the external gas at low ionization (i. e. 0°, 0+, 
S+ and N+ lines) follows the classical, Wilson's type, 
expansion law Vi,xp(km s~"'^)=15.0(±2) xR" (valid for 
1.875"<Rzt,pc<2.1G"), whereas the internal, higher ioniza- 
tion species (0++, Ar++, He+ and S++) move faster, but 
their velocity quickly decreases outwards according to the 
(rough) relation yexp(km s-^)=-43(±5) xR" -|-110(±10), 
valid for 1.725"<R2^pc<l-875". The contact is represented 



F. Sabbadin et al.: Planetary Nebulae: theory vs. practice 



15 



by H+, characterized by V;xp(H+)=28.5(±1.0) km s^^ and 
R^«pc(H+) = 1.875(±0.050) arcsec." 

To be noticed: in the internal regions of BD+30°3639 
being T4xp7^AxR", the radial depth of the zvpc changes 
across the nebula and is given by Eq. (11), where T4xp(A) 
refers to the cspl. 

The general results for the zvpc of BD+30°3639 at 
PA=10° are shown in Fig. 7, containing the line fluxes 
across the whole nebula (in arbitrary scale and corrected 
for c(H/3)=0.45; top panel), the physical conditions (T'o([N 
II]) and Ne{[S II]) at the corresponding flux peak) and the 
ionization structure of oxygen (through Eq. (32) and as- 
suming ^^^^0 for I(O",6300)>0, I(O+,7320)>0 and 

I(O++,5007)=0; middle panel) and the relative N^e]^^ 
profile ( through Eqs. (17) to (30); bottom panel). As ex- 
pected, Ha (for ro(H+)=8200 K) and 0°-hO+-hO++ (for 
re(O°)=8000 K, re(O+)=8200 K and Te(O++)=8500 K) 
provide the density distribution across the whole ionized 
nebula, whereas He+ (for ro(He+)=8500 K) and Ar++ 
(for Te(Ar"'"+)=8500 K) peak in the internal regions at 
higher ionization. 

When combined with Afc([S II]), the sharp matter pro- 
file in Fig. 7, bottom panel, has 7Ve(top)~ 15 000 cm~'^ 
(for e;=l) and the ionized mass of the nebula results to 
be Mio„~0.05(±0.01) Mq. 

Although a deep analysis of BD-|-30°3639 (whose 
multi-color and opaque movies are shown in the WEB 
page http:/ /vifeb. pd.astro.it/sabbadin) is beyond the 
aims of the paper, we wish to mention that: 

- the faint halo embedding the bright nebula (see the 
bottom panel of Fig. 7) could be a "recombining" 
halo, since it essentially retains the kinematical prop- 
erties observed in the external layers of the main neb- 
ula (vice versa, the kinematics of an AGB-halo is in- 
dependent on the main nebula, simply reflecting the 
ejection law in the pre-superwind phase; for details, 
see Sabbadin et al. 2005). If this is the case, the past 
UV flux of the exciting star was much higher than the 
present value; less than 200 years ago, the star abruptly 
cooled and/or dimmed and recombination became the 
dominant process in the outermost nebular regions, no 
longer reached by stellar radiation; 

- the combination of peculiar expansion velocity field + 
sharpness of the gas radial profile joins BD-|-30°3639 
with NGC 40, a second very-low excitation FN pow- 
ered by a low temperature and high luminosity star of 
late Wolf-Rayet spectral type (Sabbadin et al. 2000), 
and identifies a small sub-class of PNe whose internal 
layers are likely shaped by wind interaction. 

5. Discussion 

Given a regularly expanding mass of ionized gas, the zero- 
velocity-pixel-column (zvpc) of high-resolution spectral 

In other words, BD-|~30°3639 is a rare case of PN exhibiting 
the "U" -shaped expansion profile mentioned in Sect. 1. 



images represents the tangentially moving matter at the 
systemic velocity and identifies a sharp radial portion of 
nebula in the plane of the sky, whose characteristics (den- 
sity, temperature and ionic profiles) are obtainable from 
line fluxes. 

The one-dimensional distribution provided by the zvpc 
at each PA - combined with the kinematical properties 
given by the central-star-pixel-line (cspl) - represents the 
starting point for bi-dimensional tomography, i. e. recov- 
ery of the ionized gas structure within the entire slice of 
nebula intercepted by the spectrograph slit. 

In his turn, bi-dimensional tomography is the basis for 
the spatial reconstruction of the whole nebula, through a 
3-D rendering program assembling the tomographic maps 
at different PA. 

This is, in nuce, the complete logical path of our origi- 
nal methodology - valid for all types of expanding nebulae 
- overcoming any misleading camouflage due to projec- 
tion and providing the true distribution of the kinematics, 
physical conditions and ionic and chemical abundances at 
unprecedented accuracy. 

In this paper we have extracted the precious informa- 
tion stored in the zvpc of high-resolution spectral images 
of PNe, focusing on the gas density profile given by differ- 
ent atomic species. The results for hydrogen, helium and 
heavier elements can be synthesized as follows: 

(I) hydrogen (the most abundant element) presents a 
mono-electron atomic structure with a low ionization 
potential (IP=12.598 eV); at first glance, the recom- 
bination lines of H I (whose emissivity is a weak func- 
tion of the electron temperature) are excellent diag- 
nostics of the matter distribution, since they are emit- 
ted across the whole ionized nebula. However, their ef- 
fectiveness and reliability are reduced by three severe 
blurring sources - thermal motions, fine-structure and 
expansion velocity gradient within the PN - implying 
the accurate spectral image-reconstruction illustrated 
in Sect. 4.1; 

(II) the same general considerations apply to the recom- 
bination lines of both ionization stages of helium; 

(III) for heavier elements, we are forced to adopt coUi- 
sionally excited emissions of suitable ionic sequences, 
because of the faintness of recombination lines, large 
stratification of the radiation across the nebula and 
incomplete ionic coverage of the echellograms. We se- 
lect O'^+O^-l-Q"'""'" for low-to-medium ionization layers 
(IP range from to 55 eV) and Ar+++Ar'^++Ar''+ for 
the medium-to-high ionization ones (IP range from 28 
to 75 eV). The main weaknesses of this choice (and 
the corresponding care) are: (1) coUisional excitation 
rates of forbidden lines present a strong dependence 
on the electron temperature (in many cases, even a 
rough, ad hoc profile is preferable to the assumption 

Please, note the affinity between the zvpc given by high- 
resolution slit spectroscopy and the rest frame of imaging 
Fabry-Perot interferometry. 
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Fig. 7. Real proto-PNc and low-cxcitation "born-again" PNe. General properties in the zvpc of BD+30°3639 at 
PA=10° (apparent minor axis). Top panel: flux distribution of the main nebular emissions across the whole nebula 
(in arbitrary scale and corrected for c(H/9)=0.45). Middle panel: Te([N II]) (circles) and iVe([S II]) (squares) at the 
corresponding line peak and ionization structure of oxygen (from Eq. (32)). Bottom panel: relative N^ey^ profile from 
H+, He+, 0°+0++0++ and Ar++ (through Eqs. (17) to (30)). As in the case of NGC 6572 (Fig. 3), the zvpc provides 
two independent radial profiles of BD+30°3639 at opposite directions. The central star position is ar r=0.0 arcsec and 
the slit orientation is shown in the top panel. Same symbols as Fig. 3. 



Te=constant across the whole nebula) and (2) ioniza- 
tion species higher than Ar'*+ dominate in the inter- 
nal regions of high-excitation PNe powered by a hot, 
T, >100000 K, central star (A3425 A of [Ne V] , IP 
range 97.1 to 126.2 eV, could be a better diagnostic 
for the innermost layers). 

The illustrative examples presented in the previous sec- 
tions show that the "best observational" density profile is 
always a compromise among the results given by different 
elements, taking into account the general characteristics 



of each target (evolutionary phase, optical thickness - or 
thinness - and excitation degree of the nebula and tem- 
perature and luminosity of the central star). 

All this, combined with a high-quality (i. e. RR>5 and 
SS>5) spectroscopic material and a careful reduction pro- 
cedure - and thanks to the moderate dependence of A'e on 
the observational parameters (see Eqs. (17) to (30)) -, for 
the first time provides the accurate recovery (to within 10 
to 20%) of the radial matter distribution into an expand- 
ing, ionized mass of gas. 
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Our zvpc analysis 

(A) is based on the pixel-to-pixel measurement of velocity 
and flux in all nebular emissions at high, medium and 
low excitation present in long-slit echellograms cover- 
ing a wide spectral range. This strongly differs from 
(i) the usual observational setup for extended objects, 
inserting an interference filter to isolate a single echelle 
order and (ii) the normally adopted reduction proce- 
dure, aiming either at the kinematics or at average line 
fluxes; 

(B) is limited to the optical region (i.e. the ionized 
gas); for a complete understanding of the PN- 
phenomenology, it should be extended to other spec- 
tral domains, like near-IR (for the outermost envelope, 
rich in molecular and photo-dissociation emissions) 
and UV (for very-high ionization species expected in 
the innermost layers, where the star-nebula interaction 
occurs) ; 

(C) applies to other classes of expanding nebulae (e. 
g. Nova and Supernova Remnants, shells around 
Population I Wolf-Rayet stars, nebulae ejected by 
Symbiotic Stars, bubbles surrounding early spectral- 
type Main Sequence stars etc.) covered at adequate 
"relative" spatial (SS) and spectral (RR) resolutions 
(SS=R/Ar, R=apparent radius, Ar=seeing-|-guiding; 
RR=14xp/AV, AV=instrumental spectral resolution). 

The zvpc methodology developed in this paper - added to 

tomography and 3-D recovery (sec Ragazzoni ct al. 2001, 
Turatto et al. 2002, Benetti et al. 2003, Sabbadin et al. 
2004, 2005) - overcomes the long-standing de-projection 

hurdle and opens new unusual ways to deepen many facets 
of the nebular research, in particular: 

(I) to solve the topical problem of PNe - i.e. forma- 
tion, shape and shaping - by comparing the obser- 
vational results for a representative sample of targets 
at different evolutionary phases with the expectations 
coming from current theoretical evolutionary models 
(Schonbcrncr ct al. 1997, Steffcn et al. 1998, Marigo ct 
al. 2001), detailed hydro- dynamical simulations (Icke 
et al. 1992, Frank 1994, Mellema 1997, Perinotto et al. 
2004b) and updated photo-ionization codes (Ferland 
et al. 1998, Ercolano et al. 2003). At present, the most 
striking discrepancies between model-PNc (Perinotto 
et al. 2004b) and true-PNe (Sabbadin et al. 2005 and 
references therein) concern the radial distribution of 

(a) kinematics: Wilson's law is reversed in model-PNe 
(i. e. the internal layers of the bright main shell expand 
faster than the medium- to-external ones), 

(b) matter: in model-PNe the sharp density peak of 
the main shell abruptly falls inwards, whereas in true- 
PNe the whole density profile and the inwards decline 
are smooth. 

We aim at overcoming these discrepancies (both quali- 
tatively ascribable to an overestimation of the shaping 
effects due to wind interaction) by adjusting the the- 
oretical evolutive parameters until model-PNe agree 
with true-PNe; 



(II) to mine the nature and phenomenology of "exotic" 
structures at large and small scales, like haloes, caps, 
ansae, jets, blobs etc.. An emblematic case is repre- 
sented by the FLIERs (fast, low-ionization emitting 
regions) of PNe; the acronym, introduced by Balick 
(1978), identifies a series of external, N overabun- 
dant knots at low excitation and peculiar kinemat- 
ics. Although the origin and evolution of FLIERs have 
been interpreted in many fanciful ways (instability 
zones at the ionization front, condensations accelerated 
by photo-evaporation or wind-interaction, high-speed 
stellar or nebular bullets, magnetically focused knots 
etc.; Perinotto et al. 2004a and references therein), 
our high-resolution spectral survey shows that most 
FLIERs are normal condensations within the outer en- 
velope (same kinematics and chemical abundances), 
whose spectral characteristics are fully explainable in 
terms of mere photo-ionization by the central star;-'^^ 

(III) to disentangle the source (or sources) of general, 
still unsolved problems involving emission nebulae, like 
the heavy-element abundance and electron tempera- 
ture discrepancies obtained from optical recombina- 
tion and coUisionally excited lines. The proposed so- 
lutions (Peimbert et al. 2004, Liu et al. 2004, Ercolano 
et al. 2004, Sharpee et al. 2004 and references therein) 
include (a) strong temperature and density fluctua- 
tions within a chemically homogeneous nebula, (b) 
high-density, low-temperature knots overabundant in 
heavy elements embedded in a hydrogen-rich medium 
at lower-density and higher-temperature and (c) other 
physical processes (they are listed in footnote 1) com- 
peting with single electron capture (i. e. recombina- 
tion) in exciting high-level lines. 

All these (and other) hypotheses can be deeply tested 
by comparing the detailed spatio- kinematics of recom- 
bination and forbidden lines in different ions provided 
by the zvpc analysis, tomography and 3-D recovery 
applied to very-deep echellograms. 



6. Conclusions 

The presence of an expanding nebular gas is the typ- 
ical signature of the instability phases in stellar evolu- 
tion, characterized by a large, prolungated mass-loss rate 
(PNe, Symbiotic Star nebulae, shells around Wolf-Rayet 
and LBV (luminous blue variable) stars), or even an ex- 
plosive event (Nova and Supernova remnants). 

Mass-loss of evolved stars occupies a strategic ground 
between stellar and interstellar physics: it raises funda- 
mental astrophysical problems (e. g. origin and structure 
of winds, formation and evolution of dust, synthesis of 
complex molecules) , plays a decisive role in the final stages 
of stellar evolution and is crucial for galactic enrichment 
in light and heavy elements. 



For a deep analysis of the FLIERs, ansae and caps of NGC 
7009, see Sabbadin et al. (2004). 
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So far, image de-projection represented an unsour- 
montable obstacle for the accurate recovery of the spa- 
tial structure in real nebulae, thus impeding any detailed 
comparative analysis with theoretical models. Such a frus- 
trating situation is now overcome by tomography and 3-D 
recovery developed at the Astronomical Observatory of 
Padua (Italy). 

The radial density reconstruction illustrated in this pa- 
per represents the umpteenth proof that properly secured, 
reduced and analysed spectra at high-resolution constitute 
a peerless tool for deeping the spatio-kinematics, physical 
conditions, ionic structure and evolution of all classes of 
expanding nebulae, thus bridging the present, wide and 
pernicious gap between model- and true-nebulae. 
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